Signalling noise enhances chemotactic drift of E. coli 
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Noise in transduction of chemotactic stimuli to the flagellar motor of E. coli will affect the random 
run-and-tumble motion of the cell and the ability to perform chemotaxis. Here we use numerical sim- 
ulations to show that an intermediate level of noise in the slow methylation dynamics enhances drift 
while not compromising localisation near concentration peaks. A minimal model shows how such 
an optimal noise level arises from the interplay of noise and the dependence of the motor response 
on the network output. Our results suggest that cells can exploit noise to improve chemotactic 
performance. 



The motion of Escherichia coli consists of a series of 
"runs" , where the cell swims in a roughly constant di- 
rection, and "tumbles" , during which the cell randomly 
reorients [T] . In a spatially- varying environment the bac- 
terium performs chemotaxis by biasing this random mo- 
tion in the direction of favourable conditions. A well- 
studied signalling cascade [2] detects environmental lig- 
and stimuli and regulates the stochastic switching of the 
rotary flagellar motors between counter-clockwise (run) 
and clockwise (tumble) rotation pQ. However, the bio- 
chemical reactions making up this signalling pathway are 
inherently random events, and the noise introduced to 
the signal in this way will affect the swimming behaviour 
and ability of a cell to respond to gradient stimuli. Here 
we show how signalling noise can be beneficial for chemo- 
tactic performance. 

Chemoattractant stimuli are detected by the bind- 
ing of ligands to membrane receptor complexes, which 
suppresses the activity of the receptor-associated kinase 
CheA. Consequently, the phosphorylation level of the re- 
sponse regulator CheY, which in its phosphorylated form 
(CheYp) binds to the flagellar motor and promotes tum- 
bling, decreases leading to longer runs. Conversely, a re- 
duction in ligand binding leads to an increase in CheYp 
and shorter runs. Hence, the random walk motion of 
the bacterium is biased in the direction of increasing 
chemoattractant concentration. Crucially, the chemo- 
taxis network also includes a negative feedback from 
CheA activity to receptor methylation. This ensures 
adaptation of the network response to constant stimuli, 
enabling the network to detect temporal derivatives of 
the observed ligand signal [3] and allowing sensitivity to 
a wide range of ligand concentrations [4] . 

Receptor methylation and demethylation reactions are 
a significant source of noise in the signalling network [5] , 
because (i) the timescale of methylation (~ 10s [6]) is 
much longer than the other timescales in the network 
(ligand-receptor binding and receptor activity changes: 
~ 1ms; CheY phosphorylation and motor switching: 
< Is), meaning that the downstream network cannot in- 
tegrate out slow methylation fluctuations; and (ii) methy- 
lation occurs at a small number of sites on each recep- 
tor catalysed by a small number of the enzymes CheR 



and CheB, meaning that small- number fluctuations in 
the overall methylation level can be significant. Impor- 
tantly, the output of the noisy signalling network also af- 
fects the ligand signal experienced by the cell via modula- 
tion of the tumbling dynamics and hence the swimming 
trajectory. This highly non-linear feedback potentially 
means that noise may significantly affect the chemotac- 
tic response. 

Previous studies have shown that, in the absence of 
chemoattractant gradients, noise in (de)methylation re- 
actions can lead to a power-law distribution of run du- 
rations [SJ [7] . The resulting super-diffusive motion may 
enhance search efficiency compared to Brownian motion 
[8] . It has also recently been shown that slow fluctua- 
tions in the methylation dynamics can allow for enhance- 
ment of drift in linear gradients [HI H] at the expense of 
the ability of cells to localise in regions of high ligand con- 
centration [8 . However, the mechanism by which noise 
enhances drift remains unclear. 

In this paper we study the effects of receptor 
(de)methylation noise on chemotactic performance. We 
show that below a threshold noise level the steady-state 
performance of cells in a sinusoidal ligand field does not 
improve as noise is reduced. We also find that an optimal 
noise level, comparable to this noise threshold, maximises 
the drift velocity in exponential gradients. Analytic ap- 
proximations of a minimal model reveal that in the rel- 
evant regime where motor switching is fast compared to 
adaptation, drift enhancement results from the interplay 
of noise with the response of the motor switching rate to 
the output of the signalling network. 

To study the effects of signalling noise on the chemo- 
tactic behaviour of E. coli we performed simulations of 
bacterial populations using a scheme coupling swimming 
and signalling [10 . Signalling dynamics are simulated ac- 
cording to a recently-developed model j6j [11] describing 
the CheA activity, receptor methylation level, and phos- 
phorylated CheYp level. The stochastic switching of the 
two-state motor with CheYp-dependent switching rates, 
and the run-and-tumble dynamics of the bacterium in a 
three-dimensional environment are also simulated. Noise 
was introduced into the deterministic model considered 
in [TU] by adding a Gaussian white noise source to the 
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FIG. 1: (a) Mean steady-state ligand level {[L]) for a popula- 
tion of cells in a sinusoidal ligand profile with [Lo] = 800/iM 
and A — 0.25, for different noise strengths 7 -1 . Dot- 
ted lines show the maximal and minimal concentrations at 
[L] — [Lq](1 ± A). For short wavelengths A < 150/^m chemo- 
tactic cells perform worse than a simple random walk, for 
which {[L]) — [Lo] (dashed line). Increasing the noise re- 
duces the ability of cells to localise in regions of high ligand 
concentration, (b) Transient drift velocity in an exponential 
ligand profile, with [Lo] = as the noise strength 7 -1 

is varied. In shallow gradients xq > 500^im drift is fastest at 
intermediate noise strengths. 



dynamics of methylation. The strength of this noise is 
proportional to a parameter 7" 1 which is used to control 
the impact of noise without changing the response time 
of the network (full details of the model can be found in 
the Supplementary Information). 

We first considered the ability of cells to localise in 
the vicinity of ligand concentration peaks by studying 
the steady-state distribution of cells in a sinusoidally- 
varying ligand profile, [L(x)] = \Lq[{\ + A cos 2nx/X). 
Figure [iji shows the average ligand level ([£]) = 
^c^Us SiGceiis[^( x «)] experienced by the simulated pop- 
ulation as the gradient length-scale A and noise strength 
7 _1 are varied. Non-chemotactic cells performing an un- 
biased random walk would have ([L]) = [Lo]. We first 
consider cells with deterministic signalling, 7 -1 — > 0; 
the trajectories of such cells are still noisy since motor 
switching and reorientation angles during tumbles remain 
stochastic. For gradient length scales shorter than the 
typical run length of unstimulated cells, A < 5/im, cells 
are unable to react to the extremely rapidly changing 
ligand level and effectively perform an unbiased random 
walk. For long length scales A > 200fj,m cells are able 
to reliably localise in the vicinity of ligand concentration 
maxima. However, in an intermediate range of length 
scales we find that ([L]) < [Lo]. Here, methylation is too 
slow to keep the CheYp level in the sensitive range of 
the motor as cells run in directions of increasing ligand, 
and as a result cells repeatedly overshoot ligand concen- 
tration peaks, an extreme example of the 'volcano effect' 
which has been observed previously |12j . 

The shape of ([L]) as a function of A is largely un- 
changed as the methylation noise strength 7 -1 is in- 
creased; however, the difference between ([L]) and [Lo] 
at a given A is decreased. Noise moves cells away from 
concentration maxima when A is large, but also out of 
minima when A < 200/xm. However, since chemotaxis re- 



mains detrimental in rapidly-varying profiles, it appears 
unlikely that signalling is optimised for this type of en- 
vironment. We therefore focus on the regime of long 
wavelengths. Importantly, here we find that there is an 
effective noise threshold around 7 -1 ~ 0.01; provided sig- 
nalling noise is kept below this level, it is not detrimental 
for localisation near ligand maxima, which is instead lim- 
ited by the intrinsically random run-and-tumble motion 
of the cell. 

Noise reduction in biochemical signalling is typically 
energetically costly, requiring for example increased pro- 
tein production or more rapid turn-over. Our observation 
of a noise threshold suggests that there is no benefit, at 
least in terms of localisation performance, to reducing 
noise below this level. Moreover, it raises the possibility 
that this noise could somehow be exploited by the cell. 
Bacterial chemotaxis has two conflicting goals [T3]: lo- 
calisation in regions of high chemoattractant, and rapid 
drift in favourable directions. The observed chemotactic 
response has been interpreted as a compromise between 
these two objectives [T3l [14] . It is therefore important 
to also consider the effects of noise on the transient drift 
rate of cells. 

We therefore investigated transient chemotactic 
drift in an exponential ligand gradient, [L(x)] = 
[io] exp(a;/xo). Figure [TJd shows the drift velocity, es- 
timated from the linear regime of the mean x-position 
(x(t)) = N~l b X« G coiis X i( t ) as a function of elapsed time, 
for a population initially located at x = 0. We consid- 
ered only shallow gradients, xo > 500/Um, for which (x(t)} 
reaches a stable constant drift regime before saturation of 
ligand binding. Based on experimentally-observed ramp 
responses [5], cells are expected to be sensitive to gra- 
dients with xo beyond 10000/im. Interestingly, we see 
that the drift velocity has a non-monotonic dependence 
on methylation noise strength: drift in a cell with noisy 
methylation dynamics can be faster than in cells with 
no methylation noise, and there is a steepness-dependent 
noise strength that maximises the drift velocity. This ef- 
fect is also observed taking other measures of drift perfor- 
mance such as the maximal drift velocity or the mean po- 
sition of the population after some fixed time. Our results 
are consistent with previous reports that cells with sig- 
nalling noise drift more rapidly in linear gradients [51 [§] • 
We note that the optimal noise strength j^ 1 m 0.01—0.02 
gives rise to a coefficient of variation in the CheYp level 
of around 0.2, comparable to the variability required [7] 
to reproduce the experimentally-observed power law run- 
time distribution [5] . Importantly, we also see that at the 
optimal noise level neither the drift velocity in steep gra- 
dients nor the steady-state localisation performance of 
cells is significantly compromised. 

To understand the origin of this noise-induced drift 
enhancement in shallow gradients we studied a minimal 
model of chemotactic drift, shown in Fig. [2^, (see the Sup- 
plementary Information for derivation). We consider mo- 



tion in one spatial dimension; cells can be in two states, 
'+' or ' — ', which correspond to motion with velocity 
v = ±1 in directions of increasing or decreasing ligand 
concentration respectively. We assume that directional 
changes are instantaneous and hence no tumbling state 
is considered. Furthermore we assume that the CheYp 
level simply tracks CheA activity, such that the internal 
dynamics of the signalling pathway can be represented 
by a single state variable, 5a(t), that represents the devi- 
ation of the pathway activity from its adapted value (but 
with the direction of action of the gradient reversed) and 
evolves according to 5a = vr — 5a + oT(t). Here r rep- 
resents the stimulus from the ligand gradient, a is the 
strength of methylation noise, and T(t) is a Gaussian 
white noise process with unit variance. We take r to be 
a constant since in the full chemotaxis network the ef- 
fective stimulus strength goes as d x \og[L{x)} ~ x^ 1 for 
an exponential gradient. Since the noise-enhancement 
of drift is observed in shallow gradients we focus on the 
regime of small r < 1, a range of stimuli comparable in 
the rescaled parameter space to the gradients for which 
drift enhancement is observed in the full model. Finally, 
the switching propensity between '+' and '— ' states is 
given by u>(5a) — k(1 — 5a) for 5a < 1 and w(5a) = 
otherwise. This form for uj approximates the highly non- 
linear dependence of the motor switching propensity on 
the CheYp level [T5] with k setting the typical rate of 
reorientations relative to adaptation, and is chosen for 
computational simplicity. However, similar results are 
observed with w a smooth decreasing function of 5a. Fol- 
lowing |16j . evolution equations can be written for the 
joint probability p±(5a; t) of a cell to be in the '+' or '— ' 
state with internal variable 5a, 

dtp±-d Sa [{SaTr)p±} = —d 2 Sa p±Tu{5a){p + -p-) (1) 

where the first term on the right-hand side of Eq.[l]is due 
to the noise in the dynamics of 5a. Thus in each state, 
cells effectively diffuse in a potential V±(5a) ~ (5a =p r) 2 
with diffusion constant a 2 /2. The net population drift 
velocity in this minimal model is given simply by (J{t)) = 

f-oo [p+ (<*o; *) - p- ( Sa > *)] dSa - 

In the absence of noise, a = 0, Eq. [T] can be solved 
to find the steady-state drift. While the full expression 
is uninformative (see Supplementary Information), for 
small r <C 1 it can be approximated as ( J) « r/(l + 2k) 
(see Fig. J2Jd). We see that increasing k decreases (J), 
a result that also holds when a > 0; intuitively, if the 
typical run duration is shorter, less information can be 
extracted about the current direction during a single run, 
and hence the reliability with which the tumbling propen- 
sity can be regulated is reduced. 

An exact solution to Eq. [TJwhen a > is less straight- 
forward. Since in the full chemotaxis network the typical 
switching rate is fast compared to the adaptation time, 
with an effective value of ft w 10, we focus in the re- 
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FIG. 2: (a) Schematic depiction of the minimal model. 
Cells switch between the two states '+' and ' — ' with rate 
cj(8a). The internal coordinate 8a diffuses within the poten- 
tial V±(Sa). (b) Steady-state drift (J) as switching rate k is 
varied, for r = 0.25 and a — 0. Points show simulation re- 
sults; dashed line shows (J) = r/(l + 2k). (c) Steady-state 
drift (J) as the noise strength a is varied for different stimulus 
strengths r, with k = 10. For small r, (J) is maximised for 
a > 0. Lines show the approximate results for k — > oo, Eq. [3] 
(d) Probability densities p±(8a) (solid lines) when r — 0.5, 
with k = 10. For small a, p±(Sa) are similar and peaked 
around Sa — 0. For larger a, the different forces experienced 
when tj(8a) is small lead to differences between the distri- 
butions p±(Sa). Dashed lines show the effective potentials 
V cB (8a < 1) and V±(8a > 1). 

mainder of the paper on the regime k ^> 1 for which 
both numerical solutions and an analytic limit solution 
are possible. Figure [2]; shows the numerically-evaluated 
(J) at k = 10 as the noise a is varied for different stimulus 
strengths r. We can see that the minimal model qualita- 
tively reproduces the results of the full model for shallow 
gradients, with (J) showing a maximum at a > 1. 

To understand the origin of this noise-induced maxi- 
mum it is useful to consider the limit of rapid switching, 
k — > oo. In the region 5a < 1, cells rapidly exchange 
between '+' and ' — ' states, such that p+{5a) — p-(8a). 
In this region, therefore, cells spend equal time in each 
of the two states while moving in the effective mixed 
potential V cS {5a) = [V+(5a) + V-{5a)}/2 - {5a 2 +r 2 ). 
However, any cell which crosses into the region 5a > 1, 
where u>(Sa) = 0, will experience only the potential asso- 
ciated with its current state, V±(5a), until returning to 
the boundary 5a = 1. The average net drift can be cal- 
culated in terms of the mean time spent in each region, 
5a < 1 and 5a > 1 in either the '+' and '— ' states as 

(j\ = it i= (2) 

\ ' rp8a>l _|_ 2^^ a <l -)- ^<5 a >l ^ ' 

where T Sa<1 is the typical time spent in the region 5a < 1 
accounting for both the '+' and ' — ' states and we have 
used the fact that these cells have no bias in their direc- 
tion. Evaluating the typical times spent diffusing in the 
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appropriate potentials we find 

exp(cr~ 2 ) erfc(<7 _1 ) 



(J) 



V 71 
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which can be seen to have a maximum at a 2 (Fig. [2j;, 
see Supplementary Information for an exact expression 
and full derivation). 

When the noise level is small, a < 1, cells occupy the 
minimum of V e g(Sa) at 5a — and p + (5a) « p_(5a) (see 
Fig. [2ji), such that the net drift (J) w 0. As the noise 
is increased, so too is the rate at which cells reach the 
transition point 5a = 1, beyond which ui(Sa) — 0. Im- 
portantly, the offset between the minima of V+(5a) and 
V_(5a) means that cells in the region 5a > 1 in the '— ' 
state will experience a stronger force in the — 5a-direction 
than cells in the '+' state, |9,5 a V+| < |<9,5 a y_|. Hence the 
mean time spent diffusing in the region 5a > 1 before 
returning to the boundary at 5a = 1 will be longer in the 
'+' state than in the '-' state, T s + a>1 > T s _ a>1 . This is 
the origin of the drift enhancement by noise: as a is in- 
creased beyond unity, cells spend an increasing amount 
of time in the region 5a > 1, and so the magnitude of 
this effect and hence the drift (J) increase. For even 
larger values of the noise, the difference in the amount 
of time spent in the '+' and '— ' states decreases again, 
because now diffusion dominates over the difference in 
forces. Hence (J) decreases again for a ^> 1. 

While Fig. [2£ shows qualitative agreement between 
Eq.[3]and the numerical results for n = 10, (J) calculated 
for k — y oo underestimates the drift when k is finite (see 
also Supplementary Information). With a finite switch- 
ing rate, p±(5a < 1) need not be precisely identical; in- 
deed, an effective positive feedback acts on differences 
between p±(5a < 1) due to the variation of uj(5a) with 
5a < 1. As 5a becomes larger and u>(5a) decreases, cells 
will spend more time in the potential V±(5a). Since cells 
in the '+' state tend to drift towards larger values of 5a 
than cells in the '— ' state, cells will typically remain in 
the '+' state for longer, allowing for further drift and am- 
plifying the differences between p + (5a) and p- (5a). This 
means that (i) cells with 5a < 1 also contribute positively 
to (J); and (ii) more cells enter the region 5a > 1 in the 
'+' than the ' — ' state, further increasing (J). 

The maximal drift at the optimal noise strength for 
large k 3> 1 remains less than r, the drift which is 
achieved for small n — >• and a — >• (see Fig. [2J}), 
suggesting the E. coli could instead enhance chemotactic 
performance simply by reducing the switching rate and 
signalling noise. However, chemotactic performance also 
depends on steady-state localization. This motivated us 
to study the effects of varying the rate of motor switch- 
ing in the full model. As shown in Fig. [3j reducing the 
switching rate increases the transient drift velocity at low 
noise levels, but this is accompanied by a decrease in 
steady-state performance. This data emphasises that the 
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FIG. 3: Effects of varying the noise strength 7 -1 and typical 
rate of motor switching in the full model on (a) steady-state 
(A = 500/im), and (b) transient (xo = 1000/im) chemotactic 
performance. At intermediate switching rates, the optimal 
drift velocity can be achieved with a finite noise level without 
compromising steady-state performance. The estimated wild- 
type parameters are indicted with +. 



choice of chemotactic network parameters entails a trade- 
off between steady-state and transient performance. The 
switching rate cannot be significantly decreased without 
severely compromising steady-state performance. But at 
the switching rates observed in real cells, a moderate level 
of noise improves drift performance compared to a system 
without noise, without harming steady-state localisation 
and additionally reducing the cost of signalling. 

The principal result of our manuscript is an expla- 
nation of the mechanism by which noise can enhance 
chemotactic drift. This mechanism is unlike the noise- 
driven motion of Brownian ratchets or motors [17] in 
that drift reflects differences in the mean occupancy of 
internal states of the system, rather than motion in an 
underlying spatial potential. The effect described here 
can potentially enhance the bias of any two-state system 
in which the rate of switching between these states is de- 
pendent on internal variables of the system and is fast 
compared to the relaxation in either state. The essence 
of the mechanism is that the fast switching obscures the 
difference between the two states, unless noise is suffi- 
ciently strong that the system transiently passes into a 
regime of slow switching. Then the asymmetry between 
the two states is felt, and a bias is induced in the steady- 
state occupancy of the states. This effect also differs 
from stochastic resonance |18j since it is the enhancement 
by noise of a steady-state response to a constant input, 
rather than the enhancement of a dynamic response to a 
time-varying external signal. 

Additionally, we have seen that a finite level of sig- 
nalling noise can enhance transient chemotactic drift in 
shallow gradients while not significantly compromising 
drift in steep gradients or the ability to localise near con- 
centration peaks. More generally, our results suggest that 
cells may be able to exploit intracellular signalling noise 
to improve behavioural responses. 
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manuscript. This work is part of the research program of 
the "Stichting voor Fundamenteel Onderzoek der Materie 
(FOM)", which is financially supported by the "Ned- 
erlandse organisatie voor Wetenschappelijk Onderzoek 
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SUPPLEMENTARY INFORMATION FOR 
"SIGNALLING NOISE ENHANCES 
CHEMOTACTIC DRIFT OF E. COLI" 

DESCRIPTION OF THE FULL CHEMOTAXIS 
MODEL 

Here we outline the full model of the chemotaxis path- 
way, adapted from [10] . used in the simulations. Param- 
eter values used in the simulations are listed in Table 

E2 

It is assumed that the probability that a given receptor 
complex is in the active state is in quasi equilibrium and 
is given by 



1 



l + exp(iVe(m, [L])) ' 



(SI) 



where e is the free-energy difference between the ac- 
tive and inactive states and N is the number of re- 
sponding receptor dimers. The free energy e(m, [L]) = 
f m (tti) + f[L] ( [L] ) is a sum of contributions dependent on 
the receptor methylation state f m (m) = o(mo — to) and 

the ligand concentration f [L ]([L}) = -log 

where Ka and K] are the dissociation constants of lig- 
and to the active and inactive receptors respectively, 
and roughly define the range of ligand concentrations to 
which the cell can adapt. 

The dynamics governing receptor methylation and 
demethylation is given by the stochastic equation 



dm 
lit 



= k R (l — a) - k B a + r](t). 



(S2) 



Here we have introduced the Langevin noise term r)(t) 
to represent noise in receptor (de)methylation reactions, 
with (ri) = and (r)(t)ri(t')) = ~f- 1 {k R {l~a) + k B a)5{t- 
t'). The values of the rate parameters k R and k B are set 
from the experimental measurements 6 of the adapted 
kinase activity a = k R /(k R + k B ) — 0.5 (we shall 
henceforth denote values of the signalling variable in the 
adapted state with an overbar) and the slope of the feed- 
back transfer function {j^{^t)}a=a ~ ~(k R + k B ) = 
—0.03s -1 . Since the noise ri(t) is purely additive it does 
not alter the response time of m(t), which is set by k R 
and k B . The parameter 7 -1 is expected to scale with the 
number or receptor clusters as N. r ~* [5] , while the precise 
value of 7 -1 will also depend on the cell volume and non- 
diffusive transport rates of CheR and CheB enzymes. 

To propagate the noise from receptor modification, we 
include explicitly the dynamics for the fraction of CheY 
which is phosphorylated, y p , 



Parameter 


Value 


Source 


N 


6 




a 


1.7 


m 


mo 


1 


m 


K A 


3000/xM 


.10. 


Ki 


18. 2 fiM 


m 


a 


0.5 


m 


kft 


0.015s -1 


A 


kn 


0.015s -1 


m 


Vp 


0.3 


OS 


kz 


2s -1 


KB 


ky 


1.7s -1 


calculated 


b 


0.25 


® 


TO 


0.2 s 




H 


10 




P 


282250s -1 


calculated 


D e 


0.123rad 2 s -1 




V 


16.5fim.s~ 1 


union 



TABLE SI: Parameter values 

Here the first term represents the phosphorylation of 
CheY by active receptor-Che A complexes, and the sec- 
ond term accounts for dephosphorylation of CheYp by 
the phosphatase CheZ. The fraction of CheY which is 
phophorylated in adapted cells is yp — kyaj '(kyci+kz) ~ 
0.3 PS- Note that we do not introduce an additional 
source of noise from the CheYp (de) phosphorylation re- 
actions. This is because we expect such fluctuations to 
be fast compared to methylation noise, and of small am- 
plitude due to the relatively high copy numbers of CheY 
and CheZ, and hence to have a small impact on swim- 
ming behaviour. 

The tumbling bias is a sigmoidal function of y p with a 
Hill coefficient H w 10 [15] . CheYp modulates only the 
probability of the motor to switch from counter-clockwise 
to clockwise rotation, and does not affect the rate of the 
reverse transition |19j . Hence the propensities for a cell 
to switch from running to tumbling and vice versa are 
given by [TP] , 



(S4) 



dy p 
dt 



k Y a{\ - y p ) - k z y p , 



(S3) 



where r is the average duration of a tumbling event, 
which is independent of y p . The constant /3 is set 
such that the clockwise bias in the adapted state is 
b = ui r ^t(yp) I \^r^t{yp)+^>t^r) = 0.25 [HJ. Upon switch- 
ing from tumbling to running, a new orientation for the 
cell is chosen randomly and independently of the previ- 
ous run direction. Introducing correlations between the 
directions of consecutive runs, as has been observed ex- 
perimentally [1], does not significantly affect our results. 

Simulations are initialised with a population of 10000 
cells located at x = adapted to the local environment. 
Each cell initially has a random orientation and motor 
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state chosen according to the adapted tumbling bias, b. 
In a simulation step of length St, running cells move a 
distance St ■ v in the direction of their current orienta- 
tion, where the swimming speed v is taken to be con- 
stant; tumbling cells do not move. During runs, the run 
direction is also perturbed by rotational diffusion with 
diffusion constant Dg ; a random angular displacement is 
added to the orientation at each time step to account 
for this. The internal state is also updated according 
to equations ( ST]|S3 1, and the motor state switches with 
probability St ■ u; r _> f or St ■ u; t _> r , given by Eq. |S4[ as 



appropriate. 

DERIVATION OF THE MINIMAL MODEL 

Here we present the complete derivation of the minimal 
model and its parameters in terms of those of the full 
model. 

We start by considering the internal dynamics of the 
activity a. Taking the time-derivative of Eq. |Sl| and sub- 
stituting in Eq. |S2| gives 



da dm da d[L] 



da 



dt dm dt d[L] dt 



Na(l - a) 



«(&r(1 - a) - k B a + rj) + 



Kj-Ka 



d[L] 



([L]+K A )([L]+Ki) dt 



(S5) 



r 



Making the approximation Kj <C [L] <C Ka, and using 
the identity a = kji/(k R + ks), this expression simplifies 
to 



da 
~di 



Na(l - a) 



akji 1 d[L] 

— { a-a) + a V ~ W] — 



(S6) 



To model the reduced swimming dynamics in one spa- 
tial dimension we write ^4 — -i-^ffil^g^ where | e ff 
denotes an effective steepness calculated by integrat- 
ing | ^1 1 over a uniform distribution of run angles in 
three-dimensional space. For an exponential gradient 
|X(x)] = [Lq] e?q>(x/xo) we therefore have 



da 
~dl 



Na(l - a) 



ak R v 
^^{a - a) + arj =F ~ — 
a Zxq 



(S7) 



Turning to the switching dynamics, we first ne- 
glect the presence of the tumbling state since tumbling 
events are relatively short compared to runs. Then the 
rate of changing direction is simply given by u>(y p ) = 
uj T ^t(y p )/2, where the factor of 1/2 appears because 
only half of tumbling events will lead to a change of di- 
rection. Next we make a quasi-steady-state assumption 
that, since the dynamics of y p (t) is fast compared to ait), 
y p {t) = a(t)/ [a(t) + kz/ky}- Substituting into u>(y p ) and 
linearizing about the adapted value a leads to 



co(a) 



1 — (a — a) 



a{kya + kz 



(S8) 



We define the displacement variable Sa according to 
Eq.[S8]such that uj(5a) = Pfp" (1 - Sa)/2. Next we sub- 



stitute the corresponding expression for a(t) into Eq. S7 
and rescale the units of time according to t = tNak R (\ — 
a). Finally, assuming that the noise and stimulus terms 
as small deviations of the same order as Sa and retaining 



only first-order terms in Sa yields 
vHkz 



dSa 

dt 2akR(kya + kz)xo 



-So 



Hk ? 



k R (k Y a + kz) 



(S9) 

We identify the first term with vr and the rescaled ve- 
locity v = ±1, such that for a given value of r the equiv- 
alent value of xq using the parameters of Table |S1| is 
xq w 2200/im/r. Similarly, we can use the expression 
for the methylation noise strength in the full model to 
identify 



Hk, 



12(1 



kya + kz 



57 7 " 



-1/2 



(S10) 



With these parameters, the stimuli in the r < 1 regime 
of the minimal model are comparable to gradients in the 
range for which non-monotonicity in the drift velocity 
is observed in the full model. The optimal noise level 
in the minimal models is somewhat smaller than for the 
full model (a ps 2 is equivalent to 7 _1 s» 0.0012), suggest- 
ing that the linear form of uj(5a), which decreases more 
rapidly than uj r ^ t , leads to a slight overestimate of the 
effect of noise. 

It remains only to convert the mean switching rate into 



the rescaled time units, k 



2Nak R (l-a) 



10. 



EXACT SOLUTION FOR a = 

The full solution of the minimal model in the absence 
of noise is 



(J) 



Fi(tt(l+r),l+2K;-4rK) 



1 /i _ »0 i^jVli±+LL±+±l 
1 V 1 ) ,Fi(re(l+r),2re; 



r < 1 



-4r«) ■ ~ ! ( gll ) 

r > 1 



where iFi(ci, Ci\x) is a Kummer hypergeometric func- 
tion. 
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EXACT SOLUTION FOR k ->• oo 

The exact solution to the minimal model in the limit 
k — y oo can be expressed in terms of the typical time 
spent in the domains 5a < 1, and 5a > 1 in the '+' and 
'— ' states, according to 



(J) = 



T 



Sa>l 



_ 22 1 <5a<l _|_ rp8a>l 



(S12) 



To determine the typical time spent in the region 8a < 1, 
diffusing in the potential V c g(5a) — (6a 2 + r 2 )/2, before 
reaching the boundary at Sa = 1 we calculate the mean 
first-passage time from a position Sa = 1 — e, where e is a 
small displacement, to 5a = 1, using the standard result 

m 




100 



FIG. SI: The simplified Eq. 3 in the main text (full line) 
and the exact result Eq. S19 (dashed) for k — > oo. For small 
r < 0.25 the results are indistinguishable. 



T 



8a<l 



dx e 2V ^^ 2 / dy e - 2V ^l° 2 



6-S-O 



1-e 



^./^e^fl+erf 1 

G \ G 



(S13) 
(S14) 



where we have used the fact that the diffusion constant in 
the <5a-coordinate is g 2 /2. Similarly, the time spent in the 
5a > 1 region starting from a position 5a = 1 + e can be 
calculated using the potentials V + (5a) = (5a ~ r) 2 /2 + r, 



T 



8a>l 



1+e 



dx e 2V +^/° 



dy e 



-2V+(y)/c 



V G Z G 



and V- (5a) = (5a + r) 2 /2 - r, 

•l+e 



(S15) 
(S16) 



T 



8a>l 



dx e 2V ~^/° 



dy e 



-2V-(y)/a 2 



e->0 



.&( 1 ^) 2 /- 2 er fci± r 



(S17) 
(S18) 



The constant offsets to V±(5a) are included so that the 
potential landscape is continuous at 5a = 1, but have no 
effect on the final result. While each of the first-passage 
times vanishes for e — > 0, as the population of excur- 
sions into the relevant domain becomes dominated by 
extremely short trajectories which return to the bound- 
ary almost immediately, the ratios of these times remain 
well-defined. This is because the relative times spent 
in each domain are determined predominantly by long 
trajectories with macroscopic durations, rather than the 
increasing number of vanishingly-short trajectories. 
Combining the results above leads to an exact solution 



for the drift in the limit k — > oo, 



(J) = 



r erfc 



e 

l-r 



' erfc 



l+r 



2e" 



fl + erfil 



e^erfci±^ 

(S19) 

Figure SI shows the excellent agreement between the full 
result and the simplified Eq. 3 in the main text. 

Importantly, Eq. |S19| can also be derived by calculating 
the equilibrium occupancy of the potentials according to 
the Boltzmann distribution. We define the probability 
weight in the 5a < 1 region as 



W 



8a<l 



5a e - 0v °«^ 



(S20) 



with ft = 2/g 2 , and similarly for W[ a>1 . By comparison 
to Eq. |S13[ it can then readily be seen that W Sa<1 ~ 
T Sa<1 /e, and the net drift follows straightforwardly. 



[1] H.C. Berg and D.A. Brown. Nature, 239:500-504, 1972. 
[2] H.C. Berg. E. coli in Motion. Springer- Verlag, New York, 
2004. 

[3] J.E. Segall, S.M. Block, and H.C. Berg. Proc. Natl Acad. 

Set. USA, 83:8987-8991, 1986. 
[4] V. Sourjik and H.C. Berg. Proc. Natl Acad. Set. USA, 

99:123-127, 2002. 
[5] E. Korobkova, T. Emonet, J.M.G. Vilar, T.S. Shimizu, 

and P. Cluzel. Nature, 428:574-578, 2004. 
[6] T.S. Shimizu, Y. Tu, and H.C. Berg. Mol. Syst. Biol, 

6:382, 2010. 

[7] Y. Tu and G. Grinstein. Phys. Rev. Lett, 94:208101, 
2005. 

[8] F. Matthaus, M. Jagodic, and J. Dobnikar. Biophys. J., 

97:946-957, 2009. 
[9] T. Emonet and P. Cluzel. Proc. Natl Acad. Set. USA, 

105:3304-3309, 2008. M.W. Sneddon, W. Pontius, and 

T. Emonet. Proc. Natl Acad. Set. USA, 109:805-810, 

2012. 



8 



[10] L. Jiang, Q. Ouyang, and Y. Tu. PLoS Comput. Biol, 

6:el000735, 2010. 
[11] Y. Tu, T.S. Shimizu, and H.C. Berg. Proc. Natl Acad. 

Set. USA, 105:14855-14860, 2008. 
[12] D. Bray, M.D. Levin, and K. Lipkow. Curr. Biol, 17:12- 

19, 2007. 

[13] D.A. Clark and L.C. Grant. Proc. Natl Acad. Sci. USA, 

102:9150-9155, 2005. 
[14] A. Cclani and M. Vergassola. Proc. Natl Acad. Sci. USA, 

107:1391-1396, 2010. 
[15] P. Cluzel, M. Surette, and S. Leibler. Science, 287:1652- 

1655, 2000. 



[16] R. Erban and H.G. Othmer. SI AM J. Appl. Math., 
65:361-391, 2004. 

[17] P. Reimann. Phys. Rep., 361:57-265, 2002. 

[18] L. Gammaitoni, P. Hanggi, H. Jung, and F. Marchesoni. 
Rev. Mod. Phys., 70:223-287, 1998. 

[19] U. Alon, L. Camarena, M.G. Suretter, B.A. Areas, and 
Y. and Liu. EMBO J., 17:4238-4248, 1998. 

[20] C.W. Gardiner. Handbook of Stochastic Methods for 
Physics, Chemistry and the Natural Sciences, 2nd edi- 
tion. Springer- Verlag, Berlin, 1985. 



